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Abstract 

We consider systems of ODEs with the right hand side being Laurent polynomials in 
several non-commutative unknowns. In particular, these unknowns could be matrices of 
arbitrary size. An important example of such a system was proposed by M. Kontsevich. 
We prove the integrability of the Kontsevich system by finding a Lax pair, corresponding 
first integrals and commuting flows. We also provide a pre-Hamiltonian operator which 
maps gradients of integrals for the Kontsevich system to symmetries. 
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1 Introduction 



In connection with the theory of non-commutative elliptic functions, M. Kontsevich [2 J consid- 
ered the following discrete map 

U^-UVU^ 1 , V — > U^ 1 + (1.1) 

where u,v are non- commutative variables (in particular, n x n- matrices of arbitrary size). His 
numerical computer experiments have shown that this map could be integrable (see [6]). In the 
abelian case the element 

h = u + v + u' 1 + v~ l + u~ x v~ l (1.2) 

is an integral for the mapping ( 11. ip . The equation h = const defines a family of elliptic curves. 
In the non-abelian case the element h is transformed as h — > uhu^ 1 . It follows from this formula 
that trace(h k ) is a first integral of (II. ip for any natural k. 

Kontsevich also observed that ( ll.ip is a discrete symmetry of the following non-abelian ODE 
system: 

u t = uv — uv^ 1 — v^ 1 , v t = — vu + v u^ 1 + u^ 1 (1.3) 
and conjectured that (11.31) is integrable itself. 

Our paper is devoted to the system (ll.3p . It belongs to the class of systems of the form 

Ut = Pi (it, v), v t = P 2 (u,v), (1.4) 

where Pi are elements of the associative algebra A4 of all non- commutative polynomials in u, v, 
u~ l ,v~ l with constant scalar coefficients. The elements of Ai are called non-abelian Laurent 
polynomials. 

In papers [31 [1] integrable systems of type (11.41) with Pj being non-abelian polynomials in 
u, v, were considered. The existence of an infinite series of infinitesimal symmetries was taken 
as a criterion for integrability. A similar approach to integrability of evolutionary polynomial 
non-abelian PDEs was developed in [S] . As far as we know integrable systems with non-abelian 
Laurent right hand sides were not considered before. 

In this paper we find a Lax representation with a spectral parameter for system ( 11. 3p . The 
corresponding Lax L-operator generates infinitely many integrals of motion for ( 11. 3ft . They are 
integrals for the discrete map ( 11. ip as well. We also find a pre-Hamiltonian operator that maps 
gradients of first integrals to symmetries. This proves that (11.31) is integrable in the sense of 

i- 



2 Symmetries and integrals 

Let us briefly recall the definitions from [3] generalized to the Laurent case. 

Let Xi, . . . ,xn be non- commutative variables. We consider ODE systems of the form 

^ = F a (xi, x N , a^ 1 , x^ 1 ), (2.5) 
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where F a are Laurent polynomials. We denote by A4 the associative algebra of all Laurent 
polynomials. Formulas ( 12. 51) together with = — a;" 1 ^^" 1 define the corresponding 

derivation D t on A4. 

An (infinitesimal) symmetry for (12.51) is a system 
dx 

— ° = G a (x 1} ...,x N ,Xx , ...,x^ 1 ), G a eM (2.6) 

compatible with ( 12.5ft . Compatibility means that the derivations Z?t and _D T corresponding to 
(12.5 ft and ( 12. 6 j) commute. 

The existence of an infinite series of symmetries is typical for integrable non-abelian ODEs. 
The simplest non-trivial symmetry for ( 11. 3p is given by 

m t = —mm — Mf 2 + uv + (vu)" 1 + v ~ 2 + u 2 v _1 — uvu' 1 + uv~ 2 + (vuv)^ 1 + u{vuv)~ l , 

U T = fUf + fit 2 — V U — («f ) _1 — It 2 — W 2 M _1 + VUV^ 1 — VU~ 2 — (uvu)^ 1 — v{uvu)~ l . 

The following conjecture is concerned with the dimensions of vector spaces Sf. of all symmetries 
with right hand sides being polynomials in w, t> , -u -1 , t> -1 of degree < k. 

Conjecture: 

dim S4 n = 2n 2 , dim S^n+i = 2n 2 + 2n, dim 5*4 n+ 2 = 2rt 2 + 2rt + 1, dim S4„ + 3 = 2n 2 + 4n + 1. 

We have verified the conjecture for Sk, k = 0, ...,16. Notice that (11.31) is invariant with 
respect to involution 

u—tv, v — > u (2.7) 

and t — > —t. All known symmetries are either involution invariant up to r — > —t or involution 
symmetric to another symmetry. In the remainder of the paper we will denote the involution 
F\ v ++ U of any Laurent polynomial F by F. 

An element I of A4 is called a full integral for ( 12.5ft if D t (I) = 0. In the matrix case this 
means that every component of I is an integral of motion. For equation (ll.3p the element 

I = uvu^v- 1 (2.8) 

is a full integral. This element is an integral for the mapping (ll.ip as well. We were unsuccessful 
in finding more Laurent full integrals for (I1.3P (other than polynomials of I, I^ 1 )- 

In the g-case, which is a specialization of our general non-abelian situation, an additional full 
integral exists. Namely, consider the associative algebra Q generated by u, v with the relation 
uv = qvu, where q is a fixed constant. In this case the element 

J = u + qv + qu~ l + v _1 + u~ x v _1 

is a full integral for ( II. 3p . This integral is a g-deformation of the element (11.21) . 

Apart from (12.81) . equation (11.31) has so called trace first integrals. If u and v are matrices, 
these integrals are given by traces of some Laurent polynomials. For instance, the traces of h k , 
where h is given by (11.2j) and fceN, are such integrals. 

For (12. 5p with x a being non-commutative symbols we define the trace(a), a G M. as an 
equivalence class. Two elements a and h of Ai are called equivalent iff a can be obtained from 



3 



b by cyclic permutations of factors in its monomials. In other words, the traces for elements 
of A4 are defined as the corresponding elements of the quotient space A4/K, where K is the 
vector space spanned by all commutators in Ai. If a — b G K, we write a ~ b. 

An element p of Ai is called a trace integral of (12.51) . if D t (p) ~ 0. Trace integrals pi and p 2 
are called equivalent if pi — p 2 ~ 0. By definition, the degree of a trace integral is the minimal 
degree of elements from the corresponding equivalent class. 

Low degree trace integrals of (11.31) can be found straightforwardly. In particular, there are 
no integrals of degrees 1 and 3, h is the only integral of degree 2, three linearly independent 
trace integrals of degree 4 are given by h 2 ,I and I -1 . Detailed information of higher degree 
integrals can be found in Section HI 

To generate infinite sequences of trace integrals for (jl.3p the following procedure can be 
applied. 

Suppose two Laurent polynomials H,A G Ai satisfy 

D t (H) = [A, H] (2.9) 

where D t is the derivation corresponding to (11.31) . The involution (12.71) generates the pair A, H 
that is also satisfying (12.91) . It follows from (I2.9P that if is a trace integral for (jl.3p . Since H k 
also satisfy (I2.9P for any natural k the elements H k are trace integrals as well. In the matrix 
case this is equivalent to the fact that the spectrum of H is preserved under the flow (11.31) . 
Thus for any pair (Hi, Aj) satisfying (12.91) the elements 

p ik = H k , and p^=H~ k (2.10) 

are trace integrals for any natural k. 

It is easy to see that if H, A satisfy (12. 9p . then for any invertible element g G Ai the 
conjugation 

H^gHg-\ A -> gAg' 1 + g t g- 1 (2.11) 

leads to another pair satisfying (12. 9p . It is clear that conjugated pairs produce the same trace 
integrals (LTTOl . 

Apart from h given by (II .2p satisfying the relation (12 .9p with A = —u~ 1 — v we found several 
more low degree pairs (Hi,Ai) satisfying (12. 9p . 

The corresponding elements Hi, i — 1, .., 11 are given by Hi = h + a,i, where 



ai = 


vu] 


= v (S 2 /- 1 - 


1), 


a 7 


= SI + a 2 + a4, 


a 2 = 


[v , u~ x v _1 ] 


= u-\ I - 


1), 


a 8 


= SI + a 2 + [v,u~ 1 v 2 }, 


a 3 = 


[v , uv^ 1 ] 


= u (5 3 /- 1 - 


1), 


a 9 


= SI + a 4 + [u" 1 , u~ x v 


a\ = 


[w -1 , v^u] 


= v-\S I - 


1), 


aio 


= S 2 I +[v-\uv], 


a 5 = 


CLi -\- O4, 






an 


= I + [u, vu^ 1 } 


a 6 = 


a 2 + a 3 , 











Here S stands for a cyclic shift of factors in a monomial, i.e. S(abc...z) = bc.za, and I is given 
by (12. 8p . Sets of Hi, Hi that are conjugate to each other are {h, H 5 , H 6 }, {Hi, H 3 , H 2 , H4}, 
{H 7 }, {H$, H 9 , H w , Hn} and the involuted versions of these 4 groups. In addition to h, A 
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three of the pairs H h A { representing these groups up to involution and conjugation are: 

Hi — u + u~ x + V^ 1 + U~^V~^ + M _1 f U, Ai = u — V + v _1 + u^v ~ 

i?7 = M + f + U~ l V~ l U + fit - ^t? — ^ + ■U -1 !' -1 + W-U _1 f A 7 = V" 1 — V, 

H\\ = u + -u -1 + v^ 1 + -u _1 f _1 + uvu^ 1 + mm -1 !' -1 , An = —u^ 1 . 



3 The pre-Hamiltonian operator 

In the previous section we found several infinite sequences of trace integrals for ( 11. 31) . Here we 
construct the sequences of corresponding symmetries via a so called pre-Hamiltonian operator. 
This operator is defined in terms of left and right multiplication operators. 

For any Laurent polynomial a G Ai we denote by L a and R a the operators of left and right 
multiplications on Ai: 

L a (x)=ax, R a (x)—xa. 

It is clear that L a b = L a L^ and R ab = RbR a - It follows from the associativity of the algebra Ai 
that R a Lb = L b R a . 

The algebra of all left and right multiplication operators is generated by L x ., L -i, R x , R x -i- 

i i 

We denote this associative algebra by O and call it the algebra of local operators. 

For any element a = a(x) G Ai we define an 1 x N- matrix a* with entries being elements 
of O by the following identity: 

-^-a(x + e<k)| e=0 = a*(<k) . (3.12) 

For example, for h from (11.21) we have 

h* = ( 1 — L u -iR u -i — L u -\R u -iR v -i , 1 — L v -iR v -i — L u -iL v -iR v -i ). 
It is easy to see that 

A(a)=a»(F), (3.13) 

where D t is the derivation associated with (12.51) and F = (F 1; ...,F N ) T is the right hand side of 
(12~5D. 

For any vector a = (ai, oat) t , a, G Ai we define the Frechet derivative operator a* as 
the N x iV-matrix with rows (ai)*, (a^)*. 

For any two vectors p = (p±, ...,pn) t , q = (qi, ...,gjv) T , Pi,qi G Ai we put 

(p,q) = p 1 q 1 H VPnQn- 

Let a(x) G .M. Then grad (a) is the vector uniquely defined by: 

-^a(x + e(k)| e=0 ~ (<k, grad (a(x))) . 

We will denote by grad Xl (a), grad XN (a) the components of the vector grad (a). It is easy to 
see that if a ~ b, then grad (a) = grad (b). This means that grad (a) is well-defined for trace 
integrals. 
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For example, for the function h given by (11.21) we have 

grad h — (1 - u~ 2 - u~ 1 v~ 1 u~ 1 , 1 - v~ 2 - v~ 1 u~ 1 v~ 1 ) T '. 

It follows from the definition of an infinitesimal symmetry ( 12. 6 j) with symmetry generator 
G = (Gi, .., Gat) t for a system (I2.5P that D t D T x. = D T D t x, i.e. 

AG = £> T F = F*(G) (3.14) 

(from (13.13P ) is a linearized equation satisfied by G, where F* is the Frechet derivative of the 
right hand side of (1231) (cf. [4]). 

An N x N matrix V with entries from O is called a pre-Hamiltonian operator for equation 
(1231) if 

A(P) = F,7> + PF:. (3.15) 
Here the adjoint operation * on Mat/v(C) is uniquely defined by the formula 

(p,Q(q)> ~ (Q*(p),q>, (3-16) 

where Q e Mat N (0), p», $ G Al 

Relation (I3.15P can be rewritten in the form 

(D t -F*)V = V{D t + Fi). (3.17) 

It can be shown (cf. [1]) that for any trace integral a of (12.51) the vector b = grad(a) satisfies 
the equation D t (h) + F*(b) = 0. Applying both sides of (I3.17P to b, we get that any pre- 
Hamiltonian operator maps gradients of integrals for (12. 5p to symmetries. 

Proposition. The following operator (cf. [3]) 

j RuRu L U L U L U L V -\- L U R V L V R U -\- R U R V \ ^8^) 

V L U R V L V L U L V R U R V R U L V L V R V R V J 

is a pre-Hamiltonian operator for equation ( II .3p . The proof consists of the straightforward 
verification of relation ( I3.15P using 

Ry R v —i L u ~\~ L u L v —iR v —i ~\~ L v —iR v —i 

—L v — L v L u -iR u -i — L u -iR u ~i ~Ru + Ry,- 1 

L v — L v -i ~Rv ~ L u -iR u -iR v — L u -iR u -i 

Ru + L v -iR v -iR u + L v -iR v -i —L u + L u -i 

computed from (I3.12p . (13.161) . 

Applying operator (I3.18P to grad h, we get (up to a factor of 2) the right hand side of 
system (II .3p . In this sense h/2 plays the role of a Hamiltonian for ( II. 3p . However, the bracket 

{a, b} =< grad a, V grad b > 

defined on traces of Laurent polynomials does not satisfy the Jacobi identity for V from (I3.18P 
and a true Hamiltonian structure for ( II. 3p is yet unknown. 
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4 The Lax pair 



Consider 



L 



( 



v 



v 1 + u \v + v l u l +u l + l 

■r 1 + \u V + V^U^ 1 + U^ 1 + y 



) 



A 



( 



v 1 — v + u Xv 
v" 1 u 



) 



(4.19) 



where A is a (scalar) spectral parameter. Then the relation 



AL = [A,L] 



is equivalent to the Kontsevich system (11.31) . Although the pairs Hi,Ai described in Section [2] 
do also satisfy an equation similar to (12. 9p . they are not Lax pairs for (II .3p since in this case 
(12. 9 p follows from (ll.3p but not vice versa. 

The Lax pair can be replaced by any equivalent one obtained through a conjugation (12. lip , 
where g is an arbitrary invertible Laurent 2 x 2-matrix. Other equivalence transformations are 
L — > -Pi(L), A — > A + P 2 (L), where Pi are polynomials with constant A-dependent coefficients, 
and arbitrary transformation A — > /(A). 

As usual, the traces tr ~L m = (L m )n + (L m )22 generate trace integrals of motion. In partic- 
ular, tr L yields v _1 + u + v^u^ 1 + u~ l + v, which is equivalent to h from (II. 2p . In contrast 
to (12. 9p each power of L gives us several trace integrals since tr L m is a polynomial in A, A -1 
with all coefficients being trace integrals. We verified that all trace integrals of degree < 12 for 
(11.31) are generated in such a way. 

Table [1] shows all integrals of degree d generated from tr L m , m < 14 that are not generated 
from tr L l , i < m. The following statements assume that all trace integrals have been reduced 
modulo lower degree trace integrals. 

Each integral is represented by a *, o or • and is located in one diagonal of table [TJ Each 
diagonal starts in a table entry which shows a number k indicating that the integral is I k . For 
the single * - diagonal is k = 0. o - diagonals have k < and start in row 3|/c|, i.e. I k result 
from tr L 3 ' fc ' and • - diagonals have k > and start in row 4fc, i.e. I k result from tr L 4fc . 

The different powers of A in tr ~L m have the following first integrals as coefficients. * - 
integrals are the coefficients of A , o - integrals of degree d resulting from L m are the coefficients 
of \( d - 2m )/ 2 . • _ integrals of degree d resulting from L m are the coefficients of \^ 2m ~ d )/ i . 

All -k - integrals are invariant under involution (12. 7p . All other first integrals come in pairs, 
one o - and one • - integral, both are involution symmetric to each other and therefore in the 
same column. 

Applying the involution ( 12. 7p to the Lax pair (14.191) . we get a dual one: 



We tried to find a 2 x 2 Laurent matrix P such that P L = const L P (under change of A) 
but we did not succeed and we doubt such a P exists. Therefore we think that L and L are not 
related by standard algebraic symmetries. L, A and L, A are different Lax pair representations 
but each one gives all trace first integrals, only for different degrees. For example, see table 1, 




(4.20) 
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tr L 3 generates J -1 as first integral and tr L generates / as first integral. L does also generate 
/ as first integral, but from tr L 4 not from tr L 3 . And that is the case with each first integral. 
Using L instead of L does generate the same table only with o and • - entries interchanged. 
Replacing A — >■ /(A) in L (I4.19p does not change the table. 

Applying the pre-Hamiltonian operator ( 13.1 8p to * - integrals gives involution invariant 
symmetries of the same degree. Applying (I3.18P to an involution symmetric pair of first integrals 
gives 2 symmetries that are one degree higher than the first integral and that are also involution 
symmetric to each other. 

Traces of Hj (I2.10p have been verified to be linear combinations of integrals of table [TJ For 
example, 

• first integrals tr h 2m+k , n — 0, 1, .., k — 0,1 are of degree 2(2m + k) and require only up 
to L 2m+fc to be derived, 

• first integrals tr(h + ai) 2m+k , m = 0, 1, .., k — 0,1 are also of degree 2(2m + k) (although 
h + ai is of degree 3) but require up to L 4m+fc to be derived, 

• first integrals tr(h + a-j) m and tr(h + au) m , m — 0,1, .. are of degree 4m and require up 
to L 3m to be derived. 



All first integrals generated by the Lax pair operator L(= L^^) are linear combinations 
of first integrals generated from L and vice versa. 

If a given first integral F of degree 2(2m + k), m = 0, 1, k = 0, 1 is to be expressed as a 
linear combination of first integrals computed from L and L then it requires at least L 2m+fc or 
L m+ and at most L 3m+fc and L m+ . If F is to be expressed as a linear combination of first 
integrals computed from L alone then it requires at most L 4m+fc . 

The table of first integrals with its straightforward extension appears to be complete because 
all symmetries of degree up to 16 have been verified to be generated from integrals of this table 
and the pre-Hamiltonian operator (I3.18p . 



5 Summary 

In this paper we provide two Lax pair representations L, A in ( I4.19P and L, A in (I4.20p for the 
non-abelian Kontsevich ODE system (II .3p . We also give a pre-Hamiltonian operator V in ( I3.18P 
that converts gradients of first integrals into Lie symmetries. By computing independently all 
symmetries with Laurent polynomial generators up to degree 16 we show that both Lax pairs 
generate all trace first integrals at least up to this degree. For any n we propose a conjecture 
on the dimension of the vector space of symmetries of order < n. 

Outstanding work contains an algebraic description of all A-operators for the hierarchy and 
a recursion operator for the model. 

Another open problem is to find a Laurent Lax pair for the original discrete Kontsevich 
map. We have checked numerically for 2x2 and 3x3 matrices u, v that the characteristic 
polynomial of L(A) does not change under the discrete map (11.11) . In this case we can apply 
the discrete map in components, it is invertible. But this is not the case on the level of Laurent 
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Table 1: # of new trace first integrals of degree <i generated from tr L m not generated from 
tr L l , i < m. 

polynomials since we cannot apply the map to L given in (I4.19p . It is an outstanding problem 
to find a Laurent matrix P such that the discrete map is applicable to L = PLP~ X . 
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